Impurity spin dynamics in 2D antiferromagnets and superconductors 
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We discuss the universal theory of localized impurities in the paramagnetic state of 2D antiferromagnets where 
the spin gap is assumed to be significantly smaller than a typical exchange energy. We study the impurity spin 
susceptibility near the host quantum transition from a gapped paramagnet to a Neel state, and we compute 
the impurity-induced damping of the spin-1 mode of the gapped antiferromagnet. Under suitable conditions our 
results apply also to d-wave superconductors. 



Doped antiferromagnets (AF) have been the 
subject of intense studies in the context of the 
cuprate high-temperature superconductors and 
other layered transition metal compounds. We 
present a quantum theory of a particular class of 
doped AF where it is possible to neglect the cou- 
pling between the spin and charge degrees of free- 
dom and consider a theory of the spin excitations 
alone. Such a theory will apply to (i) quasi-2D 
'spin gap' insulators like SrCu203 or NaV^Os in 
which a small fraction of the magnetic ions (Cu 
or V) are replaced by non-magnetic ions like Zn 
or Li and to (ii) high-temperature superconduc- 
tors like YBa2Cu307 in which a small fraction of 
Cu has been replaced by non-magnetic Zn or Li. 
In the first case the spin gap A is significantly 
smaller than the charge gap justifying a theory 
of the spin excitations alone. In the second sit- 
uation the effect of the fermionic quasiparticles 
in the superconducting state can be shown 0] to 
be weak due to the linearly vanishing density of 
states of the Fermi level. 

The effect of a (magnetic or non-magnetic) 
impurity can be probed by measuring the uni- 
form spin susceptibility, which takes the form 
X = {g^B) 2 (Axb + Ximp) where A is the total 
area of the AF, \b is the bulk response per unit 
area, and Ximp is the additional impurity contri- 
bution. In the paramagnetic ground state of the 
host each impurity induces a distortion of the host 
spin arrangement with a net magnetic moment 5 
associated with the impurity. The distortion is 
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plies that the impurity susceptibility follows 



Xi 



imp 



5(5 + 1) 
3k B T 



as T -> 0. 
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For a non-magnetic impurity in a spin-1/2 system 
we have 5 = 1/2; for a general impurity eq. (0) 
can be used as definition of 5. 

The basis of our investigations is a boundary 
quantum field theory which describes a bulk AF 
together with arbitrary localized deformations. 
We focus on the vicinity of a quantum transi- 
tion from a paramagnet to a magnetically ordered 
Neel state: Then the spin gap in the paramag- 
netic state is small compared to a typical nearest- 
neighbor exchange, A <C J, which is the situa- 
tion realized in many compounds. The field the- 
ory has been discussed in Ref. jjj; it consists of 
d + 1-dimensional (f> A theory for the bulk order- 
ing transition and a coupling to a local quantum 
impurity spin. The renormalization-group (RG) 
analysis shows that both the bulk and the bound- 
ary couplings are marginal for d = 3 and flow to 
fixed-point values for d < 3. This implies that 
the coupling between the bulk and impurity exci- 
tations becomes universal, and the spin dynamics 
in the vicinity of the impurity is completely de- 
termined by bulk parameters, the gap A and the 
velocity of spin excitations c. Based on the RG 
results one obtains a number of universal prop- 
erties in an expansion in e = 3 — d; we mention 
here the behavior of the uniform susceptibility at 
the bulk critical point. The system shows the 
Curie response of an irrational spin as T — > 0, 
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Ximp = C\/ (ksT), where C\ is a universal number 
independent of microscopic details. The e expan- 
sion result for C\ is 
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More detailed dynamic information can be ob- 
tained by a self-consistent diagrammatic method. 
The paramagnetic phase of the bulk is assumed 
to be dimerized, its spin-1 excitations can be de- 
scribed using triplet bosons t\^ a . The impurity is 
represented by an additional spin S a at site 0, 
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where J is the host exchange constant, et the 
energy of the spin-1 mode in the bulk, K the 
coupling constant to the impurity spin, and N s 
the number of lattice sites. The impurity spin is 
represented by auxiliary fermions /, the impurity 
dynamics is contained in the fermion self-energy 
which arises from the scattering off the t bosons. 
We employ a self-consistent non-crossing approx- 
imation (NCA) to calculate this self-energy; this 
approach follows from a saddle-point principle af- 
ter generalizing the spin symmetry to SU(iV) and 
taking the limit N — > oo. The NCA equations 
can be solved in the scaling limit; the value of the 
coupling K drops out of all results for physical 
observables provided that A <g; J - we obtain the 
same universal behavior as predicted by the RG. 
In fact, the results for susceptibility and impu- 
rity spin correlations agree with the one-loop RG 
result 0. 

The diagrammatic approach can be easily ap- 
plied to a system with a finite density of impuri- 
ties nimp. The important observation is that the 
impact of the impurities is determined by a sin- 
gle energy scale Y = n imp (fhc) d / A d_1 . The AF 
in the absence of impurities shows a pole in the 
dynamic susceptibility xq( w ) a t the AF wavevec- 
tor Q. Our main concern is the fate of this col- 
lective peak upon the introduction of impurities. 
Scaling arguments predict that the susceptibility 
takes the form 
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Figure 1. Left: Universal lineshape Im<!> for 
T/A = 0.125. Right: The same, but convoluted 
with a Gaussian corresponding to the experimen- 
tal resolution of || , together with the data points 
of Ref. ||. We have used A = 43 meV - this 
small shift may be attributed to perturbations 
being irrelevant in the RG sense. 



where $ is a universal function, and A denotes the 
quasiparticle weight. In the absence of impurities 
we have $(u7,0) = 1/(1 - (u; + i0+) 2 ). The self- 
energy of the spin-1 bosons caused by the scat- 
tering at randomly distributed impurities is cal- 
culated using a self-consistent Born approxima- 
tion. The equations for the Green's functions can 
be entirely written in terms of scaling functions 
with arguments Hu)/A and T/A, consistent with 
the scaling prediction (^) . A numerical result for 
$ is shown in Fig. [j]. The quasiparticle pole is 
broadened to an asymmetric line, with a tail at 
high frequencies. Our theory can be applied to a 
recent experiment Q where the inpurity-induced 
broadening of the spin-1 'resonance peak' at en- 
ergy A = 40 meV in YBa2Cu307 has been ob- 
served. This experiment has ni mp = 0.005, T = 5 
meV, T/A = 0.125. The half-width of the line 
is approximately Y, and this is in excellent ac- 
cord with the measured linewidth of 4.25 meV, 
see Fig. [j]. More tests of the predictions of our 
theory should be possible in the future. 
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